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We study the size fluctuations of a local denaturation zone in a DNA molecule in the presence 
of proteins that selectively bind to single-stranded DNA, based on a (2 + l)-dimensional master 
equation. By tuning the physical parameters we can drive the system from undisturbed bubble 
fluctuations to full, binding protein-induced denaturation. We determine the effective free energy 
landscape of the DNA-bubble and explore its relaxation modes. 



PACS numbers: 87.15.-v, 82.37.-j, 87.14.Gg 

Under physiological conditions the Watson-Crick 
double-helix is the thermodynamically stable configura- 
tion of a DNA molecule. This stability is effected by the 
specific Watson-Crick H-bonding, whose key-lock princi- 
ple guarantees the high level of fidelity during replication 
and transcription; and by the stronger base-stacking be- 
tween neighboring base-pairs causing hydrophobic inter- 
actions between the planar aromatic bases yj, |2j • 

The initial breaking of the stacking interactions in an 
unperturbed DNA molecule is associated with an acti- 
vation barrier o$ ~ 10 -3 - -5 0, 0, 0. Once this bar- 
rier is overcome, the free energy of breaking an addi- 
tional base-pair is of the order of 1 to 2k B T 0, 0, ef- 
fecting local, single-stranded DNA denaturation zones of 
a few tens of broken base-pairs \J\. By thermal acti- 
vation, these DNA-bubbles fluctuate in size by a zipper 
motion at the two forks where the bubble meets the in- 
tact double-strand [2j . By fluorescence correlation meth- 
ods this DNA-breathing can be explored on the single 
molecule level, revealing multistate (un)zipping kinetics 
with a typical single (un)zipping step time scale around 
50/j,s and a typical bubble lifetime of the order of 10 ms 
0. 

A long standing puzzle had been why the presence of 
selectively single-stranded DNA binding proteins (SSBs) 
does not lead to full DNA-denaturation, as SSB bind- 
ing is thermodynamically favorable • Detailed single- 
molecule optical tweezers studies, by overstretching of the 
DNA molecule to bring the effective temperature close to 
the melting temperature T rn in the presence of bacterio- 
phage T4 gene 32 SSBs, showed quantitative ly lV\\ that 
there exists a kinetic block for SSB-binding |l2j: below 
T m the bubble lifetime is shorter than the typical SSB 



binding time, counteracting helix-destabilization. 

In what follows, we develop a dynamical model to 
quantify the coupled dynamics between a fluctuating 
DNA-bubble and SSBs that attempt to bind to it, in 
terms of the system parameters (temperature, external 
force, SSB binding rate and strength, SSB size). We 
demonstrate that the presence of SSBs leads to enhanced 
bubble lifetime. Effectively, the bubble free energy is low- 
ered, and even SSB-induced denaturation can occur. 

The system we have in mind (Fig. ^) resembles the 




Figure 1: Clamped DNA bubble in a region of size M, im- 
mersed in a bath of SSBs. SSBs do not bind across zippers. 



DNA-construct from Ref. 0, in which a homopolymer 
bubble region is clamped at both ends. The one-bubble 
approximation used here is generally valid below T m due 
to (To ^ 1- Fig. ^ also illustrates that the typical binding 
size A of an SSB (~ 10 bases), is of the same order as 
the bubble size. It is therefore necessary to consider the 
statistical weight from SSBs explicitly, instead of defining 
an effective chemical potential |T^ . 

To quantify the system, we define the probability dis- 
tribution P{m,n,t) to find a bubble of size m, with n 
bound SSBs at time t. Its time evolution is governed by 
the (2+l)-dimensional master equation 



dP(m,n,t)/dt = t + (m- l,n)P(m- l,n,t) + t~(m + l,n)P(m + l,n,t) 

— (t + (m, n) + t~ (m, n)) P(to, n, t) 

+r + (m, n — l)P(m, n — 1, t) + r~ (to, n + l)P(m, n + 1, t) 

— (i- + (to, n) + r~ (to, n)J P(m, n, t), (1) 
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owing to the discrete nature of the problem. In Eq. , the transfer rates t ± describe changes in the bubbles size 
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to, and changes in the number n of bound SSBs. The 
boundary conditions for the are reflecting at n = and 
"max = 2 [to/ A], the maximum number of SSBs that can 
bind to a bubble of m denatured base-pairs. Similarly, 
we impose a reflecting boundary condition at to = and 
to = M, the maximum bubble size. Moreover, we have 
to consider that a bubble cannot zip close if its size is a 
multiple of the SSB-size, m = kX (k £ N) and the number 
of bound SSBs is n max or n max — 1, i.e., at least one of 
the arches of the bubble is fully occupied. To define the 
transfer rates based on the statistical weight 3f(m, n) to 
find a state (m, n), we assume detailed balance [lall6j : 



t + (m — 1, n)3f(m — 1, n) 
r + (m, n — l)3f(m, n — 1) 



t (to, n)iF(TO, n) (2) 
r~ (m, n)5°(m, n). (3) 



The statistical weight 3f(m,n) = J^to) J° SSB (m,n) 
has two contributions. The bubble part according to the 
Poland-Scheraga model for DNA-melting is 0, 03 



(4) 
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J^-(to) = cr li m (l + my c , TO > 1, 

with the bubble initiation factor <jq , the weight u 
for breaking a base-pair, and the loop closure factor (1 + 
m)~ c for creating a polymer loop of size to, with the offset 
by f due to persistence length corrections 0, 11(1. We 
choose typical values, <7o = f0~ 3 and c = f .76 |l8l l2fij|. 
Eq. Q is completed by 3?~(Q) = 1 . In comparison to the 
continuum description of DNA-breathing in absence of 
SSBs by a Fokker-Planck equation involving the gradient 
of the free energy [2lJ , we note that the discrete approach 
(|T|) explicitly includes the constant activation barrier ao, 
compare Fig. |21 

The contribution from the SSBs has the form [rHHsf 



3f SSB (m,n) = K n tl(m,n), 



(5) 



with the binding strength k = coK e<1 , involving the SSB- 
concentration Co and equilibrium binding constant K eq = 
vq exp(/3|i?ssB|), where is the typical SSB volume and 
-E-SSB its binding energy. The weight fl(m,n) counts all 
possible ways of putting n SSBs onto the two arches of 
the bubble, including potential gaps between SSBs due to 
the fact that their size A is larger than a base. Explicitly, 



fi(m, n) = oj(to, n')uj(m, n — n') 



(G) 



with n, (n — n') < ri max /2, and the combinatorial term 

/ m — (A — l)n 



uj[m, n) 



(7) 



By detailed balance and taking the rate t + for breaking 
a base-pair to be proportional to the activation factor u 
(including loop closure effects), and assuming that the 
rate r - for SSB unbinding is proportional to the number 
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Figure 2: Bubble free energy in absence of SSBs as function of 
bubble size m ( — : c = 1.76; : c = 0). We chose ao = f0~ 3 . 



n of bound SSBs, we arrive at the expression for the 
transfer coefficients, i.e., the rates for the bubble size 

t~ (to, n) = kfl(m — 1 , n)/Q(m, n) (8a) 
t + (TO, n) = ku([l + m]/[2 + to]) c , to > f (8b) 

with the bubble initiation rate t + (0, 0) = 2~ c kaou, and 



r (to, n) = n"/k 



(9a) 



r + (to, n) = fkn(n + I )0(to, n + 1 )/Q(m, n) (9b) 

for the SSB number transfer rates |22| |. Here, we in- 
troduced the dimensionless ratio 7 = q/k of the SSB 
unbinding rate q and the base pair zipping rate k. 

To solve the master equation we introduce an 
eigenmode expansion of the form 



P(m,n,t) = ^c p Qp(m,n)exp ( - t/r^j, 



(10) 



in which the coefficients c p of a given eigenmode p are 
determined via the initial conditions. The correspond- 
ing eigenvalue equation for the bubble size-SSB number 
eigenfuction Q p determines the mode relaxation times t p , 
and can be solved numerically, or, in some limits, analyt- 
ically 0| . From the bubble-size autocorrelation function 



A{t) = (m(t)TO(O)) 



p#0 



-t/r p 



(11) 



a typical quantity determined in experiment, we obtain 
the relaxation time spectrum {t p } with the correspond- 
ing amplitudes A p = (^ mn TnQ p (m,n)\ . The slowest 



mode T r , 



Ti determines the characteristic equilibra- 



tion time |23f. These measures can be used to quantify 
the bubble-SSB system for different cases: 
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Figure 3: Relaxation time spectrum for <to = 10 and M = 
40. For it = 0.999, the longest relaxation time dominates. 

(i.) In absence of SSBs, Eq. reduces to a (1 + 1)- 
dimensional master equation. If we neglect the loop clo- 
sure factor, we can obtain an analytic result using orthog- 
onal polynomials, from which we infer the inequalities for 
the eigenvalues t p 



k 



1 



,1/2 



<t m < ■■■ <n<k 1 M 

(12) 

The equal signs hold in the limit M — » oo. The charac- 
teristic relaxation time r re i ax tends to very large values 
(and diverges for M — •> oo) in the limit u — ► 1, i.e., on 
approaching T m . In Fig. [31 we plot the corresponding 
relaxation time spectrum for three temperatures. Well 
below T m , we see the multistate relaxation observed ex- 
perimentally For temperatures closer to T m , the slow- 
est eigenmode becomes increasingly dominant (w 2-state 
behavior). To estimate the rate constant k for zipping we 
find from Fig.[3]that at u — 0.6 T re i ax ~ 20/fc; comparing 
to the experimental result T rc i ax ~ 10ms, we consistently 
extract k ~ 50/iS. 

(ii.) The fast binding limit corresponds to 7 3> 1, 
i.e., SSB (un)binding being much faster than the bubble 
zipping. By adiabatic elimination [l6| . the master equa- 
tion Q reduces to an equation in the bubble size m as 
in case (i.) but where the bubble free energy landscape 
[3^ — — ln^d, with = Yin ^'(' m : n )i is dressed by 
the ongoing SSB (un)binding. This effective free energy 
is lowered in comparison to case (i.), as demonstrated 
in Fig. 21 for two cases: (a) the effective free energy is 
reduced, but the overall slope still positive; (b) the SSB- 
dynamics causes a negative slope of the effective free en- 
ergy, ultimately leading to complete denaturation of un- 
damped DNA. We note a characteristic feature of the 
effective energy in Fig. 0] namely, the finite size effects 
due to the SSB size A: the bubble has to open up to a 
minimum size m — A, before it is able to accommodate 
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Figure 4: Effective free energy in the limit 7 > 1 ( — ), and 
'free energy' for various fixed n (u — 0.6, M = 40, c = 1.76, 
A = 5). Top: k — 0.5; bottom: stronger binding, k = 1.5. 



the first SSB, etc. This also effects a nucleation barrier 
for SSB-assisted DNA-denaturation even in the case of 
negative slope of the dressed free energy for larger m. 
Numerically, we determine the critical binding strength 
K crit « 1 for the parameters of Fig.0] at which the effec- 
tive fi^(m) is almost flat. Fast binding is, e.g., realized 
for gp32 mutants at elevated effective temperature [ll| . 

(Hi.) General case. If SSB (un)binding is not suffi- 
ciently fast but occurs within the typical bubble relax- 
ation time in absence of SSBs, the full master equation 
{JU needs to be solved. By tuning the relative SSB un- 
binding rate 7 and SSB binding strength k, we can move 
from the situation with practically no binding to the fast 
binding case. In Fig. we illustrate this behavior via 
the characteristic relaxation time T rc i ax . Experimentally, 
u and k are changed by variation of temperature and 
concentration of SSBs in solution, respectively. 

To (un)zip close (open) a base-pair, the two single- 
strands making up the bubble have to be pulled closer 
towards (pushed away from) the zipper fork. This addi- 
tional effect may be included using similar arguments as 
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where particles bind to a substrate whose binding sur- 
face fluctuates in time (wetting, oxidation or binding of 
biomolecules to membranes); in particular, in the pres- 
ence of finite size effects. In addition, the coupled DNA 
bubblc/SSB binding dynamics studied here is a generic 
example of how a stochastic process (partially) can rec- 
tify another one. 

We are happy to thank Richard Karpel, Oleg 
Krichevsky and Mark Williams for helpful discussions. 
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Figure 5: Relaxation time as function of u, for k = 0.5, ao = 
10" 3 , M = 20, A = 5. 

in Ref. 0. The adjustment of pulling (pushing) prop- 
agates along the contour of the chain until a bend (in- 
flexion) is reached, a distance that scales as the gyration 
radius, i.e., ~ m v . Having in mind Rouse-type dynamics, 
this would slow down the (un) zipping rates by a factor 
m~ v [la] ■ The relevance of this effect for the rather small 
bubble sizes well below T m is not obvious, and has to be 
based on more accurate experimental or numerical data. 

At temperatures well below T m , our description will be 
valid for undamped DNA homopolymers due to <To <C 1. 
Explicit boundary and heteropolymer structure effects 
can be included by introducing the bubble position as an 
additional variable in the master equation [H]. Our 
description will also hold under moderate chain tension, 
e.g., in combination with optical tweezers setups. The 
pulling force / then gives rise to blobs of size £ = fcsT// 
in which the DNA is undisturbed j^. Conversely, in 
the case of strong pulling twist is progressively taken out 
of the DNA molecule, corresponding to a torque T and 
the twist energy 8qT, where 6q = 27r/10.35 denotes the 
twist angle per base. The statistical weight is modified 
according to u — > uexp((38 T) 

Our (2+l)-dimensional master equation approach pro- 
vides a quantitative framework for the coupled dynamics 
of the size fluctuations of DNA denaturation bubbles and 
the binding of SSBs. It explains the different regimes, 
from free bubble breathing to SSB-induced denaturation, 
that correspond to experimentally accessible in vitro as 
well as in vivo situations. In particular, we expect that 
this scheme is useful to the design of future experiments, 
and to estimate in vivo conditions based on the SSB bind- 
ing strength k, the base stacking factor u, and the ratio 
7 of SSB unbinding and base pair zipping rates. 

Given the rather general formulation in terms of a mas- 
ter equation for the probability distribution P(m,n,t), 
our approach may be useful to various other systems 
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